Tempered relaxation with clustering patterns 



Aleksander Stanislavsky'^'^, Karina Weron'^'^ 

""Institute of Radio Astronomy, Ukrainian National Academy of Sciences, 
4 Chervonopraporna St., 61002 Kharkov, Ukraine 
^Institute of Physics, Wroclaw University of Technology, 
Wybrzeze Wyspianskiego 21, 50-310 Wroclaw, Poland 



Abstract 

This work is motivated by the relaxation data for materials which exhibit 
a change of the relationship between the fractional power-law exponents 
when different relaxation peaks in their dielectric susceptibility are observed. 
Within the proposed framework we derive a frequency-domain relaxation 
function fitting the whole range of the two-power-law dielectric spectroscopy 
data with independent low- and high-frequency fractional exponents 7 and 
—a, respectively. We show that this effect results from a contribution of dif- 
ferent processes. For high frequencies it is determined by random stops and 
movement of relaxing components, and the low-frequency slope is caused by 
clustering in their temporal changes. 
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1. Introduction 

The stochastic approach based on subordination of compound random 
processes for description of the non-exponential relaxation phenomena within 
the anomalous diffusion framework shows itself as one of the most useful 
mathematical tools [1]. Its recent success allows one to get the fractional 
two-power dependency commonly observed in relaxation of complex mate- 
rials (glasses, liquid crystals, polymers, biopolymers and so on). According 
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to this approach, the well-known Havriliak-Negami (HN) relaxation appears 
due to a compound subordination by two independent random processes. 
One of them is connected with a coupling between the very large jumps in 
physical and operational times, and another accounts for the amount of time 
when a relaxing entity does not participate in motion. In fact, both of them 
are described by means of a skewed a-stable process and its inverse. Al- 
though, from one side the results help one to understand better the complex 
nature of the physical mechanisms underlying the power-law relaxation re- 
sponses, at the same time they raise a number of problems to the theory 
of relaxation. In particular, any a-stable process has even no the first mo- 
ment. In real situations there are sufficiently many factors truncating the 
distribution of a-stable processes This leads to the tempered 

a-stable processes, and they have all the moments. As it is shown in [3], 
from the subordination by the inverse tempered a-stable process one can de- 
rive the tempered diffusion equation and the relaxation function describing 
the Debye (D), Cole-Cole (CC) and Cole-Davidson (CD) types of relaxation. 
The tempered diffusion has a transient character, i. e. a crossover from sub- 
diffusion at small times to normal diffusion at long times. The transient 
subdiffusion has impact on kinetics of magnetic bright points on the Sun 
[sl and has been observed in cells and cell membranes j^, 10, 11 1. Physical 



arguments for appearance of such effects are that subdiffusion is caused by 
traps, but in a finite system there is a given maximal depth of the traps 
(maximal waiting time) truncating their power-law waiting time distribution 
in such a way that beyond the maximal waiting time the diffusive behavior 
of the complex system tends to normal. Note also that the truncation of 
waiting times demonstrates features of weak ergodicity breaking in motion 



of lipid granules [12|. The above mentioned relaxation functions are only 
partial cases of the more universal HN law (see Fig. [T]). Let us mention at 
this point that the original HN function fits 13, the fractional power- law 
dependence observed in large part of the dielectric susceptibility xi'^) data: 

^Xi^) = xiO) - xi^) oc (iu/up)"" for w < , (1) 
xi^) oc {iu /u)p)"'~^ for u ^ Up 

with < m < 1 and < 1 — < 1, the low- and high-frequency expo- 
nents, respectively. It has to be stressed the HN exponents fulfill the fol- 
lowing relation m > 1 — n. Here Up denotes the loss peak frequency, (zo;)'*' 
means (iu)^ = \u\^ exp(i A tt sgn(a;)/2) and i = ^/^. The data character- 
ized by the opposite relation m < 1 — n (the diagram part of Fig. [1] below 
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the line CC) cannot be interpreted in terms of this function, since its ori 



gins, found within the fractional Fokker-Planck equation [15| and continuous 



time random walk |16l . Il7j approaches, cannot yield this relation. The cor- 
responding time-domain description uses the notion of a response function 
f{t) = —d(f){t)/dt (negative time derivative of the relaxation function (j){t)), 
exhibiting the following power-law asymptotics 
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(t/Tp)-" for t<rp 



for t>rp. 



where Tp = l/up is the characteristic relaxation time, determined by the 
loss peak frequency. Here is useful to mention that an ideal capacitor with 
an impedance of the form l/{iujC) does not exist in nature, and as a rule, 
dielectric materials exhibit a more realistic fractional behavior l/{iujC)°' with 
a < 1 (l8|. 

The main open question, arising here, is whether it is possible to describe 
such a general relaxation picture (in the sense of arbitrary exponents m and 
1 — n) by the tempered scheme. It would be interesting and important for 
understanding of the relaxation mechanism, especially when a change of the 
relationship between the power exponents from m>l— ntomKl — nis 
observed in one dielectric system for different dielectric susceptibility peaks. 
As we see, the same material can have both the m < 1 — n and m > 1 — n 
relaxation patterns under different temperature/pressure conditions, and a 
model, in which this passage via model parameters is possible, is required. 
It would be also of a great importance to find a single analytical form of the 
relaxation function fulfilling such requirements (see Fig. [T]and Tabled]). 

The purpose of this Letter is to shed light on the character of relaxation 
processes in such materials. To reach this goal we present a new scenario 
of subordination based on finite-moment random processes that extend the 
theory of tempered relaxation [3]. The basic idea presented here also starts 
with a compound subordination of two random processes independent of each 
other. One of them is the inverse tempered a-stable process but the another 
one is the undershooting process with index < 7 < 1 [l], |l9[ . 

The Letter is organized as follows: at first, we briefly define the pri- 
mary random processes and formulate a modification of their subordination 
representation. Next, within this approach we derive the frequency-domain 
relaxation function in an analytical form. In the last section the results of 
our study are discussed as applied to relaxation phenomena. 
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Figure 1: (Color online) Relaxation diagram positioning different laws of relaxation. The 
exponents m and 1 — n focus on declination of the imaginary susceptibility x" i'^) for low 
and high frequencies. The circles are experimental points (for various materials) taken 
from the book 



Table 1: The relation between exponents m and 1 — n for some experimental data (see 
also [i3). 



Material 


m 


1 — n 


Comments 






Polyvinylide- 


0.07 


0.49 


a-peak, m < 1 — 


n 




ne fluoride 












Polyvinylide- 


0.53 


0.13 


/3-peak, m > 1 — 


n 




ne fluoride 












Polyvinylide- 


0.05 


0.12 


7-peak, m < 1 — 


n 




ne fluoride 












Nylon 610 


0.6 


0.57 


T = 483 K, m > 1 




n 


Nylon 610 


0.3 


0.51 


T = 413 K, m < 1 




n 


Nylon 610 


0.18 


0.41 


T = 373 K, m < 1 




n 


Glycerol 


0.55 


0.46 


P = 3.1 kb, m > 1 




n 


Glycerol 


0.5 


0.69 


P = 4.4 kb, m < 1 




n 
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2. Tempering and coupling in the compound subordination 



The tempered a-stable process [20|, l2l| is characterized by the foUowing 



Laplace image of its probability density function (pdf) 

/,(w) = exp(5°-(w + 5)") , (2) 

where the stability parameter < a < 1 and the tempering parameter 5 > 
are constants. If S equals to zero, the tempered a-stable process becomes 
simply a-stable. In other words, the parameter 6 provides just a truncation 
of the ordinary, long-tailed totally skewed a-stable distribution. The trun- 
cation leads to the random process having all moments finite. Formula 
describes probabilistic properties of the tempered process in terms of the in- 
ternal (operational) time. Its inverse process may be used as a subordinator. 
The pdf ga{T,t) of the subordinator depends on the real physical time and 
describes the first passage over the temporal limit t. Its Laplace transform 
reads 

ga{T,u) = --^/(u,r) = 

U OT 

^-^^-^ exp {-t[{u + 5r - 51) . (3) 

u 

The inverse tempered a-stable process accounts for motion alternating with 
stops so that the temporal intervals between them are random and with 
heavy tails in density. The main feature of the process is that it occurs only 
for small times jsj. 

If the 7-stable process is subordinated by its inverse, the compound sub- 



ordinator behaves 1-similar (proportional to time t) for long times [22| . Such 
a subordinator does not, however, lead to any fractional two-power relaxation 
dependency like the HN one. Therefore, we consider a more complex subor- 
dination, where a coupling between physical and operational times (7-stable 
process and its inverse, respectively) is directed by an independent inverse 
tempered a-stable process. The coupling leads to two random processes un- 
derestimating and overestimating the real time t: 

T-[S^{t)] < t < T^[S^{t)] for t > , (4) 

where T^(t) is the 7-stable process, S^{t) = inf{r : T^(r) > t} is its inverse, 
and T~{t) = limj.^T-_ (x) is the left-limit 7-stable process. The random 
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process X~{t) = T~[S^(t)] has finite moments of any order wliereas X^{t) = 
T^[S^{t)] has even no tfie first moment (see details in [l9j). 

We take into account tfie process witli tlie finite moments only, i.e., 
X~(t) = T~[S^(t)] with < 7 < 1. Its probability density is of the form 



p-(y,t) = '-^y^-\t-y)-\ 0<y<t. (5) 



TT 



Its moments of any order read 

Process X~(t) = T~[S^(t)] is 1-similar and evolves to infinity as the real 
time t. 

To obtain the widely observed fractional two-power relaxation pattern 
([1]) we construct a compound subordinator which combines two following 
random processes: one is the coupling process X~{t), and another is the 
inverse tempered a-stable process Qa(t). The process X~{t) directed by 
Qaif) is just the new compound subordinator, namely Wa^-yif) = X~ {Qa(t)). 
We assume that X~(t) and Qait) are independent on each other. In this 
case the anomalous diffusion is determined by the following relation 

/•oo /"OO 

p''{x,t)= / p^{x,y)p'{y,T) go,{T,t)dydT , (6) 
Jo Jo 

where p^(x, t) is the probability density of the subordinated process obtained 
from Y directed by Wa^-yit), and p^{x,y) is the probability density of the 
parent process Y. The explicit form of the latter is not important here. For 
simplicity, but without loss of generality, as the parent process can be taken 
the standard Brownian motion. The process X~{t) is the essence of time 
clustering in the compound subordination. This process is approximated by 
a simple continuous-time random walk in which each waiting time is exactly 
equal to the jump. Consequently, a walker, moving along a Brownian tra- 
jectory in presence of the subordination by X~{t) (without any government 
of Qa{t)), stops from time-to-time and overjumps through intermediate po- 
sitions in the Brownian trajectory. In another words the process X~{t) leads 
to a random partition (or clustering) of walker's trajectories. 
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3. Derivation of relaxation function 



The primary, commonly accepted, interpretation of relaxation phenom- 
ena is based on the concept of exponentially relaxing objects (for example, 



dipoles) with different relaxation rates [13||. If the relaxing objects do not 
interact with each other, then their macroscopic relaxation function is de- 
scribed by the D law (pit) = exp{—Upt). But the simplest situation, as it 
is known from many experiments, is realized very rarely. Since the relaxing 
objects can interact with their environment, their evolution to an equilibrium 
state has a complex (random walk-like) behavior which, in general, does not 
yield the D law. This interaction may be taken into account with an aid of 
the temporal subordination [g, 23[. In this case the relaxation function reads 



0(t) 



where p{ri, t) is the subordinator pdf. Particularly, when such a subordinator 
is the inverse a-stable process Sa{t), the relaxation function takes simply 
the Mittag-Leffler (or the corresponding CC) form. Note that the solutions 
of fractional differential equations are just expressed in terms of the Mittag- 
Leffler function and its various generalizations (see, for example, the excellent 
books 0,1251 )• 

They have extensive modern-day applications in the study of 
complex systems, which maintain the long-memory and nonlocal properties 
of the corresponding dynamics such as anomalous diffusion 26|] and non- 



exponential relaxation [27 
with impact phenomena 



fractional dynamics [31 



28l . |29| . |30[ . In particular, multi-particle systems 
that exhibit interactions between particles, have 
Nevertheless, the fractional differential equations 
are only a macroscopic description of the complex system evolution, and their 
microscopic dynamics should be found out from other reasons. 

If the exponentially decayed states are subordinated by the random pro- 
cess Wa^-yit), then the time-domain relaxation function is expressed in terms 
of the integral relation 



oo POO 



^0 



e ''^yp {y,T)g^{T,t)dydT 



sm 7r7 



TT 



oo /•! 



(7) 



'0 ^0 

where Up is a constant. For the experimental study the frequency-domain 
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representation of the latter function 



dm] 



is of more interest. It is well known [l3[ that the complex dielectric 
susceptibility x{^) = x'{^) ~ ix"{'^) oc v^*(w) of most relaxing substances 
shows a peak in the loss component 'x"{<^) at a characteristic frequency ujp. 
After integration of Eq. ([8]) by means of Eq. ([7]) we get 

where 0<cr = 5/a;p<cxDisa positive constant. Such a frequency-domain 
relaxation function fits the asymptotic behavior of the dielectric susceptibility 
with different low- and high-frequency power tails, namely 

Ax(cj) oc {iu /upY for w ^ , a 7^ 0, 

oc (zcj/wp)"" for w ^ . (10) 

The main feature of the susceptibility obtained above is that the indices 
m = 7 and 1 — n = a are determined by different processes. The index 
a is connected with the tempered a-stable process only, and 7 is basically 
determined by the coupling process T~[S^{t)\. Thus, by the law ([9]) we can 
describe any fractional two-power relaxation dependency with both m > 1—n 
and m < 1 — n (see Fig. [2]). This advantage allows one to get any point on 
the relaxation diagram shown in Fig. [H but here is some nuances. 

Recall that if one uses the conventional a-stable process instead of the 
tempered one, then in order to come to both fractional two-power relaxation 
dependencies, both undershooting X~{Sa(t)) and overshooting X^{Sa(t)) 
processes in the corresponding compound subordinator should be considered 
In particular, the overshooting process clearly leads to the HN law of 
relaxation, but the process itself has no finite moments. Notice that the 
ordinary a-stable process infiuences on both indices (m and 1—n) in the 
response function (see Table [2]), whereas the undershooting process acts only 
on the index m, and the overshooting process changes only the index 1 — n. 
As the indices < a,7 < 1, the value 07 < a holds true always. With 
the parameter 6 = the tempered a-stable process becomes simply a-stable, 
and the frequency-domain relaxation response ip*{uj) coincides with the form 
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Figure 2: (Color online) Real and imaginary terms of susceptibility x(w) 0, showing the 
power-law property ([1]). Dashed lines represent m — j — 0.75 and 1 — n — a — 0.6 whereas 
dash-dotted lines correspond to m = 7 = 0.6 and 1 — n = a = 0.75. The parameter a 
equals to 0.5. 



arising from the aforesaid undershooting process in the corresponding com- 
pound subordinator jsij. In this case the low-frequency asymptotics Ax{uj) 
is proportional to {iu/ujpy^ . The CC law is just obtained from Eq.(7) for 
(7 = and 7 = 1. The CC plot for Eq. ([9]) is shown in Fig. [31 

When the parameter 7 is equal to 1, the expression (j9]) simplifies to the 
form characteristic for the ordinary tempered relaxation. This is not sur- 
prising because in this case the probability density p~{y,T) becomes the 
Dirac-delta function 6{t — y). When the parameter a tends to 1, the proba- 
bility density ga{t, t) reads as 6{t — r). For 6 = a = 1 we obtain clearly the 
ordinary exponential relaxation. Moreover, the real and imaginary parts of 
the permittivity data as a function of frequency obtained for the investigated 
samples of Cdo.ggMno.oiTeiGa at various temperatures are well fitted by the 



dependence ([9]) with the parameters a ~ 0, a ~ 1 and 7 ~ 0.63 [33[. This 
dependence is also applicable for a ppr opriate glass-forming materials such 
as, for example, ones studied in [ij, |34 . 
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Figure 3: (Color online) Cole-Cole plot for Eq. ^ with the same values of parameters as 
in Fig.[2j the solid (symmetric) line corresponds to the CC susceptibility with m = 1 — n = 
0.75 and 6 = 0, the dashed (asymmetric) line relates to the tempered case for m > 1 — n 
whereas the dash-dotted line (asymmetric in another side) represents the tempered case 
for m < 1 — n. 



Table 2: Asymptotic behavior of the frequency-domain response function for fractional 



relaxation scenarios 


(see also |l]). 






Type of relax 


- Operational time 


m 


1 — n 


ation 


(subordinator) 






Mittag-Leffler 


Sa{t) 


a 


a 


(ML) 








Havriliak- 




a 




Negami 








Generalized 




a7 


a 


Mittag-Leffler 






(GML) 








Tempered 








response Qj 


X-{Qa{t)) 


7 


a 
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4. Conclusions 



This approach gives a chance to include the fractional two-power-law 
dependency, widely known from relaxation experiments, with m > 1—n and 
m < 1 — in the framework of tempered relaxation. Before, with this in 
mind to describe the diagram of Fig. [H the HN relaxation law was completed 
for m < 1 — n by the GML (Generalized Mittag-LefHer) dependence (see 
more details in jH, [ssj). Although the new form of the frequency-domain 
relaxation function is more complicated than the HN's law, the latter requires 
necessarily to account for the GML relaxation law because a considerable 
part of experimental data in relaxation studies cannot be described by any 
modified version of the HN law. The proposed attempt also permits one 
to select easily effects influencing on the power character of low- and high- 
frequency tails in the relaxation function. The appearance of the passage 
from rnKl — ntom>l— n and vice versa in the same material reflects 
structural changes to show what has a dominant character, either flip-flops of 
dipoles (relaxing objects) at short times or cluster patterns of dipole sets at 
long times. The short-time behavior can be produced by traps, and the long- 
time trend is a result of long-range interactions of dipoles via their cluster 
and collective regions. It should be stressed that in this scheme we use only 
the random processes with finite moments what is important for analysis and 
interpretation of experimental data 14 . 



Acknowledgments 

Work of K.W. was partially supported by the Ministry of Sciences and 
Higher Education project PB NN 507503539. A.S. is grateful to the Institute 
of Physics and the Hugo Steinhaus Center for pleasant hospitality during his 
visit in Wroclaw University of Technology. 



References 

[1] K. Weron, A. Jurlewicz, M. Magdziarz, A. Weron, and J. Trzmiel, 
Phys. Rev. E 81 (2010) 041123. 

[2] P. Hougaard, Biometrika 73 (1986) 387. 

[3] R.N. Mantegna, H.E. Stanley, Phys. Rev. Lett. 73 (1994) 2946. 
[4] I. Koponen, Phys. Rev. E 52 (1995) 1197. 



11 



[5] S.J. Boyarchenko, S.Z. Levendorskij, Int. J. Theor. Appl. Finance 3 
(2002) 549. 

[6] A. Stanislavsky, K. Weron, and A. Weron, Phys. Rev. E 78 (2008) 
051106. 

[7] A. Stanislavsky, K. Weron, Phys. Lett. A 373 (2009) 2520. 

[8] A.C. Cadavid, J.K. Lawrence, and A. A. Ruzmaikin, Astrophys. J. 521 
(1999) 844. 

[9] M. Platani, I. Goldberg, A.L Lamond , and J.R. Swedow, Nat. Cell 
Biol. 4 (2002) 502. 

[10] A. Wedemeier, H. Merlitz and J. Langowski, Europhys. Lett. 88 (2009) 
38004. 

[11] U. Schmidt, M. Weiss, J. Chem. Phys. 134 (2011) 165101. 

[12] J.-H. Jcon, V. Tejedor, S. Burov, E. Barkai, C. Selhuber-Unkel, K. 
Berg-S0rensen, L. Oddershede, and R. Metzler, Phys. Rev. Lett. 106 
(2011) 048103. 

[13] A.K. Jonscher, Dielectric Relaxation in Solids, Chelsea Dielectrics 
Press, London, 1983. 

[14] A.K. Jonscher, Universal Relaxation Law, Chelsea Dielectrics Press, 
London, 1996. 

[15] Y.P. Kalmykov, W.T. Coffey, D.S.F. Crothers, and S.V. Titov, Phys. 
Rev. E 70 (2004) 041103. 

[16] K.Weron, A. Jurlewicz, and M.Magdziarz, Acta Phys. Pol. B 36 (2005) 
1855. 

[17] A. Jurlewicz, K. Weron, and M. Teuerle, Phys. Rev. E 78 (2008) 
011103. 

[18] G.W. Bohannan, "Analog realization of a fractional control element 
- revisited", in Wavelength Electronics, Inc., Las Vegas, NE, USA, 
2002, pp.175-182. 



12 



[19] A. Stanislavsky, K. Weron, Phys. Rev. E 82 (2010) 051120. 

[20] A. Piryatinska, A.I. Saichev, and W.A. Woyczynski, Physica A 349 
(2005) 375. 

[21] J. Rosinski, Stochastic Proc. Appl. 117 (2007) 677. 

[22] W. Feller, Introduction to probability theory and its application. Vol. 
II, John Wiley & Sons Inc., New York, 1967. 

[23] A. Stanislavsky, Phys. Rev. E 67 (2003) 021111. 

[24] A. A. Kilbas, H.M. Srivastava, and J.J. Trujillo, Theory and Applica- 
tions of Fractional Differential Equations, Elsevier, Amsterdam, 2006. 

[25] Anomalous transport: foundations and applications, Eds. R. Klages, 
G. Radons, and I.M. Sokolov, Wiley- VCH, Weinheim, 2007. 

[26] R. Metzler, J. Klafter, Physics Reports 339 (2000) 1. 

[27] F. Metzler, W. Schick, H. G. Kihan, and T. F. Nonnenmacher, J. 
Ghem. Phys. 103 (1995) 7180. 

[28] V.V. Novikov, V.P. Privalko, Phys. Rev. E 64 (2001) 031504. 

[29] R. Hilfer, Ghem. Phys. 284, (2002) 399. 

[30] V.E. Tarasov, J. Phys.: Gondens. Matter 20 (2008) 145212; 175223. 

[31] J. A. Tenreiro Machado, Gommun. Nonlinear Sci. Numer. Simulat. 16 
(2011) 4596. 

[32] J. Trzmiel, K. Weron, A. Stanislavsky, and A. Jurlewicz, Phys. Rev. 
E 83 (2011) 051102. 

[33] J. Trzmiel, private communication. 

[34] I.-S. Park, K. Saruta, and S. Kojima, J. Phys. Soc. Jap. 67 (1998) 
4131. 

[35] A. Stanislavsky, K. Weron, and J. Trzmiel, Europhys. Lett. 91 (2010) 
40003. 



13 



